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magnetic nanostructures, the thermokinetic theory is applied to a metallic ferromagnet in contact 
with a reservoir of spin polarized conduction electrons. The spin flip relaxation of the conduction 
electrons is described thermodynamically as a chemical reaction. The diffusion equation of the 
chemical potential (or the giant magnetoresistance) and the usual Landau-Lifshitz-Gilbert (LLG) 
equation are derived from the entropy variation. The expression of the conservation laws of the 
magnetic moments, including spin dependent scattering processes, leads then to the generalized 
LLG equation with spin polarized current. The equation is applied to the measurements obtained 
^k^j ■ on single magnetic Ni nanowires. 
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An unexpected and spectacular effect due to spin polarization of conduction electrons in metallic ferromagnets, the 
, Giant Magnetoresistance (GMR) , appeared with the first transport studies on magnetic nanostructures |l| . The spin 
diffusion length of conduction electrons being of some few tens of nm, the relaxation of the conduction electron spins 
becomes observable when the magnetization can be controlled over this typical length. Some predictions about the 
inverse effect, namely the influence of spin polarized current on the dynamics of the magnetization, were also proposed. 
Berger predicted the existence of some surprising phenomena due to the action of spin polarized conduction electrons 
on domain walls or spin waves j^] in magnetic thin films. Slonczewski predicted the rotation of the magnetization 
due to polarized current in multilayered systems and Bazaliy et al. derived from microscopic considerations a 
generalized Landau-Lifshitz-Gilbert equation ||. All the above mentioned approaches are microscopic and based on 
the ballistic approximation. 

From an experimental point of view, Freitas and Berger, Hung and Berger ||, and Salhi and Berger |7j show the 
action of a high current density on domain walls in thin films. Recent experiments on nanostructured samples bring 
I ' new evidence for the interpretation in terms of the action of the spin of the conduction electrons. Tsoi et al. Q show 
, an effect of a high current density on spin waves generation in Co/Cu multilayers, Sun reported on current-driven 
1 magnetic switching in manganite ||, and Myers et al. reported an effect of current induced switching in magnetic 
multilayer device [ |10| . In a recent work, we have evidenced an effect of current-induced magnetization reversal in 
magnetic nanowires, fTl}! , fl2]| where the reversal of the magnetization is induced by a high current at an applied field 
20% smaller than the normal reversal field. The ballistic approximation is however difficult to justify in all these 
' experiments. 

A phenomenological approach based on the thermokinetic theory |lj| , jl4| of a metallic ferromagnet in contact with a 
reservoir of spin polarized conduction electrons is proposed. In contrast to the pioneering works of Johnson and Silsbce 
|l5f on spin polarized current, the spin- flip scattering is introduced here as a chemical reaction. This formulation 
allows the conservation laws of the spin polarized conduction electrons and the ferromagnetic order parameter to 
be treated (section I). This phenomenological approach models the different effects able to take place in the system 
by a set of coupled transport equations (section 11). The differential equation of the chemical potential is derived 
leading to the known formula of the GMR (section 111). For a closed ferromagnetic system without electric current, 
O ■ the common Landau-Lifshitz-Gilbert equation is derived (section IV) . The open system with both ferromagnetization 
, and spin polarized conduction electrons leads to the description of the generalized Landau-Lifshitz-Gilbert equation 
k>! ' with polarized current (section V). The consequences in terms of current-induced magnetization reversal are deduced. 
rS , The model is applied on the data obtained on Ni nanowires. 
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I. CONSERVATION LAWS 



I 



Spin-dependent transport processes in layered structures are described on the basis of the following simple picture. 
After entering in the k-th layer , the incident current (which was spin polarized along the axis described by 
the unite vector ±itk-i in the layer E fc_1 ) first aligns Jll| along the axis ±Ufc. Inside the ferromagnetic layer, the 
population of spin up (N + ) and spin down (N_) is then not conserved due to spin-flip scattering, and some of the 
down spins relax to the up direction. 
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In this picture, the states of the sub-system E fe are described by the variables 



(S k ,Al k ,N k ;.N k ,M ) 



(1) 



where S k is the entropy, Mq = M k Uk is the magnetization of the layer k without current, N± is the number of 

■-, k 

conduction electrons with spin parallel to the unit vector iiik and Mq is the time derivative of the magnetization in 
the layer k without current. 

If Mq and N± are independent variables, the conservation of the magnetic momentum reads 



dM k ^ k ,- k - fcw 

— — = M - N_)uk- 



(2) 



where \ib is the Bohr magneton and g is the Lande factor. 

In order to write the conservation laws, the spin-flip scattering mechanism is described as a chemical reaction 
transforming a spin down into a spin up along the axis ±mj.. In this context, the reaction rate 4? is introduced as 
the number of chemical events per unit of time jL7|, fisfl . Let j^* k+1 anc j jk-^k+i ^ e ^ e curren ^ f particles flowing 
from the layer S fe to the layer S fe+1 due respectively to the electrons with spin in the direction Uk and to the electrons 
with spin in the direction —Uk- The conservation of the particles is then described by 



dt 
dN* 

dt 



= a{k-l;k)Il-^ K + {l-a(k-l;k)) I K _ 
= (l-a(k-l;k)) I k -^ k + a (k-l;k)I k 



fc-i- 



k^k+1 



-l^k _ 
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where a is the spin flip probability of the alignment process |l9| ]. In the case of ballistic alignment a(k — 1; k) = 
cos 2 i\0{k — l;fc)) with the angle between Uk~i and Uk- Introducing the polarized current I p , and the normal 
current Tjv defined by 



rk-l^k 
V 

T k-l^k 
2 N 



fe-i- 



fe-i- 



j-k-l- 
jk-1- 
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the Eq. (g) can be put into the form 

dN 




7 (fc- l;k)I k - 1 ^ k +* fe 
7 (fc- l]k)I k - 1 ^ k +* fe 



(5) 



with 7(> - 1; k) = 1 - a(k - 1; k) = sin 2 (\6(k - 1; k)) 

Inserting Eq. (H) into Eq. (0) the conservation of the magnetization reads: 



dM k 
dt 



Mn 



gfiBUk-i 



-l^k 
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k^k+l 



7 (fc-l;fc)/*- 1 ^ fc + 
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The problem of the spin transfer between the polarized current and the magnetic layer is hence solved if the polarized 
current I p and the reaction rates ^ can be describe as functions of the experimentally accessible parameters, the current 
In, the electric field Eq and the kinetic coefficients. This task is typically performed by the application of the first 
and second laws of the thermodynamics. 



II. KINETIC EQUATIONS 

The system S fc is open to heat transfer, to chemical transfer, and to mechanical work due to the magnetization and 
magnetic fields. Let us define the heat and chemical power by P^, and the mechanical power by Pw- The first law of 
the thermodynamics applied to the layer T, k gives 

r>k— 1— >fe jztk^k+1 | rjext^k 
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where P^ t ^ k — —H ext -M k . Furthermore, with using the canonical definitions T k = , fi± = ^jL^ , H k = 



the energy variation is: 



dE k 
dt 



rj-ik 



,dS k 
dt 



. k dN k 
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~dt 



,dM k 



dE k dM$ 
dM k dt 



(8) 



In the present work we limit our analysis to the isothermal case, T k = T. The entropy variation of the sub-layer is 
deduced from the two last equations, after introducing the conservation laws: 
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dE k dM k 
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k dt 
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/§ = n\ + n k _ 

jjext.k = jjext , _ rpj^ c h enuca i a ffi n ity of the reaction, defined by A k = = //; 

introduced. 

The entropy being an extensive variable, the total entropy variation of the system is obtained by summation over 
the layers 1 to where the layer 1 is in contact to the left reservoir R l and the layer is in contact to the right 
reservoir R r . Letting 

A k = A k - 2g^ B (H k - 1 - H ext - k - 1 ) 



where the total chemical potential is /i§ = (A. + (JL. We have furthermore defined H 



k 1 = H k ■ u fe _i and 
— fx_ has also been 
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the total entropy variation is: 
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where the two first terms in the right hand side of the equality stand for the heat and chemical transfer from the 
reservoirs to the system E. 

The variation of entropy takes the form 



T dS = -yp.jri + pext^ 



dt 



(12) 



where Fi are generalized forces and X 1 are the conjugated generalized fluxes. The variation of entropy is composed 
by an external entropy variation P ext (t)/T and by an internal entropy variation dS mt /dt. 

By applying the second law of thermodynamics dS mt /dt > we are leading to introduce the kinetic coefficients 
l a p such that dS mt /dt = J^i-^i (X/j hj^J ■ By identification with the expression (fy]), the kinetic equations are 
obtained: 



r fe-i->fe 

N 
P 

-I k 
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dMft 
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(13) 
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The indices N and p stand respectively for the normal and polarized transport processes (see section III), the 
indices c stands for the spin-flip scattering chemical reaction and the indices M and M account for the dynamics of 
the magnetization (see section IV). 

The kinetic coefficients are state functions: Uj = lij(S k , M k , N+, TV*, M ) and the symmetrized matrix is posi- 
tive : i {Iji + hj}uj\ > 0. Furthermore, according to Onsager relations, the kinetic coefficients are symmetric or 
antisymmetric hj = Hji 

We assume in the following that cross effects between electronic transport and ferromagnetic transport are negligible, 
so that lij = if i = {N,p, c} and j = {M, M}. Note that a polarized current is directly produced by a non-uniform 
magnetization state, through the coefficient ^(k — I, k) in Eq (|l3|). 

The physical meaning of the kinetic coefficients is described in the two following sections. 



III. GIANT MAGNETORESISTANCE 



We focus in this section on the first three equations of ( |l3| ) which describes the electric transport with spin 
polarization. After performing the continuum limit, we have: 



JjV 
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(14) 



where 7 = liuidz^o ^ and the choice of symmetric coefficients L p n — Ln p and L^ c ~ L c n is motivated by the 
fact that we neglect in this work the direct effects of the magnetic field on the charge carriers |2(| . 

In the framework of the two-channel approximation pl| ] , the coupling between the two conduction bands is neglected 
(i.e. there is no cross effect between the currents J+, J_, and 



L 



Nc 
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(15) 



where <jq > is the mean conductivity of the two spin channels. The conductivity asymmetry (3 of the two channels 
is given by 



L 



Np 



Equation Eq (14) leads then to the set of equations: 
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(16) 



(17) 



with 1 > 1 and L cc > 0. 

In the stationary state, ^gf- = 0, and assuming that (3, cto an d 7 a re approximately independent of z the diffusion 
equation of the chemical affinity is deduced (see appendix): 



d 2 A 

dz 2 
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L DW 



A + k 
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where the spin diffusion length l s f is given by 

Isf = y 

the "Domain Wall" diffusion length low is given by 



00 (1 - P 2 ) 



2eL c 



(19) 



DW 
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(20) 
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and the parameter k is given by 



k s ^<T$) < 21 > 

Note that the chemical affinity A is equal to the difference of the chemical potentials of the two conduction bands 
A = /!+ — /!_. If we assume no rotation of the spin polarization axis, 7 = (which implies antiparallel magnetic 
configuration at the interface), then equation Eq ( |l8| ) is the well known diffusion equation describing the so-called 
"spin accumulation" or "spin depletion" effect responsible for the giant magnetoresistance |15| |, p3[ , Q, A 
straightforward calculation (see appendix) leads to the giant magnetoresistance of the interface: 

flGMR _ 2__i^___ l af (22) 



IV. LANDAU-LIFSHITZ-GILBERT (LLG) EQUATION 

In this section we focus on the magnetic transport equation without electric current: Jn = Jp = 0. From Eq (|^) 
the entropy variation reduces to: 

T dS = m _ pl + (# _ &ext) dMo + (_j>E \ (23) 

at dt \ qj^ I dt 



So that the application of the second law of thermodynamics yields, 

I TJ fjext\ _ T dMg i f dM ( 

[M - H ) - Imm -jT + 1-mm — 



( dE \ _ i dM , f dMp 



(24) 



where the kinetic coefficients l a p are the coefficients of the inverse matrix {l a p\ 



{a/3}' 



Note that in adiabatically closed systems, f — f^~J > H an< ^ are functions (i.e depend only of the state 

variables (S, Mq, Mq), and not on H ext ). Since the kinetic coefficients are also state functions, the first equation in 
@ shows hence that ^ depends on H ext . We are then leading to impose l MM = in order to satisfy the second 
equation in (|^) , which gives the magnetic kinetic energy J2(| . The coefficient l MM can be identified to the magnetic 
mass, and the first equation in (E4j) gives the total magnetic force F mag acting on the system: 



b = MM ~j£T ~ K H ~ H ) - l MM (25) 



equation (25) rewrites: 



fmag = JLf E _ jj eXt .M ) - T] — (26) 

dM at 



where we have identified the Gilbert friction coefficient [£7|, to ?y = Imm- 
The theorem of the kinetic momentum gives the equation of the dynamics: 

^ = r (Mo x F-A = TM x (- - ,^>) (27) 
dt \ J \ 8M dt j 

where T is the gyromagnetic ratio and the magnetic Gibbs potential j?9) is defined by V = E + Mo-H ext ■ Equation 
( f27| ) is the well known Gilbert equation ^7j, pSjl , and can be put into the following Landau-Lifshitz form. In the case 
of uniform magnetization we have Mo = M s uo, where M s is the saturation magnetization, Eq (p7[) rewrites 
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U 



?' (u x w) - h'u x (ti x w) (28) 



where V is here the gradient operator on the surface of a unite sphere. The phenomenological parameters h' and 
are linked to the gyromagnetic ratio T and the Gilbert damping coefficient 77 by the relations |5(J 

h> - Ta 

11 ~ (l+a 2 )M s 

a = ?7rM s 



V. LLG EQUATION WITH SPIN POLARIZED CURRENT 

Let us assume an interface composed by an incident current 7* of conductivity asymmetry —j3, polarized in the 
direction e p which enters in a ferromagnetic layer (F) polarized in the direction u with conductivity asymmetry (3. 



dN 1 : 



The transfer of magnetic moments is describes by the term J = — N_ . Equation (|5|) rewrites 



dt 



= i;- I* - 2jl£ = (3 eff I N (29) 



where f3 e ff = 2/3(1 — 27(1 — 7)) and the expressions of I* and P p are derived in the appendix. 

The change of the magnetic moment of the layer due to the polarized current is given by equation (^|) and (|28| ) : 

u w - g' (u x Vy) - h'u Q x (u x w) + ^ p ef f I N e p (30) 

where the first, second and third term in the right hand side are respectively the precession term (or transverse 
relaxation), the longitudinal relaxation term, and the spin transfer due to spin polarized conduction electrons. 




FIG. 1. Uniform magnetization and the magnetic field in the case of uniaxial anisotropy. 



In order to estimate the effect of the injection of spin polarized current, the equation (|30|) is applied to the case of 
monodomain ferromagnet with applied field oriented at the angle 9 from a single anisotropy axis (see Fig. ^). If the 
vector u makes an angle (p from the anisotropy axis, the Gibbs energy density can be written in the following form 

V(ipip) = KS{-cos 2 tp - 2h(cos(6)cos(tp) 

+ sin{6)sin{(f)cos{ip)) (31) 

where h = H ext /H a is the reduced applied field defined with the anisotropy field H a , K = H a M s is the anisotropy 
constant, S is the section and ip is the out-of-plane coordinate of the vector u. Due to the cylindrical geometry, tp = 0. 
Before injecting the current, the angle ipo is given by the equilibrium condition VT^ = 0. 



G 



The precessional term can be neglected in (|30|) (low frequency response and/or high damping limite pOU ), and 
M w M . 

Experiments and samples are described in Ref. jl2| , |32| , |33| . Ni nanowires are obtained by the method of elec- 
trodeposition in track etched membrane templates. A micro-contact is realized, and the magnetoresistance of a single 
nanowire is measured. The wires are about 70 nm diameter and 6000 nm length and the magnetic energy is domi- 
nated by the Zeeman energy term and the shape anisotropy (or magnetostatic term) , very close to that of an infinite 
cylinder. The anisotropy field is calculated to be naH a « 0.3 T. 

The effect of the spin-polarized current was evidenced experimentally by injecting a strong current of about 
2. 10 7 A/ cm 2 at a fixed value of the external field h — h sw — Ah smaller than the field h sw where the switching 
occurs without current. The magnetization switch occurs at the angle tp c {9). The maximum distance Ah where the 
jump of the magnetization can still be observed corresponds then to the variation of the angle Aip = ip c — ipo needed 
to shift the magnetization up to the unstable state. 

For steady states, inserting h = h sw — Ah, ip = ip c , Eq. ( ^o|) leads to 

Ah = h sw (0) ~ 2 ^(^;)--(^ C ) (32 ) 

where v is the polar vector perpendicular to u. The parameter c is defined by the relation 

Peff h~ 
c = — — — 

eKvaCn 



(33) 



where the activation volume v a of magnetization M s was estimated to be v a ss 10 22 m 3 , K « 10 5 J/m 3 p2| , and 
/3 eff w (3 w 0.3 ||, a rj 0.15 §§. We obtain c w 200A" 1 . 

All parameters in Eq. (|3^) are known if the magnetization reversal mode, which describes the irreversible jump, is 
known. In some few theoretical models of magnetization reversal fjlf , the functions H sw (9) and ip c (0) are analytical. 
In the framework of the present empirical approach, the experimental data are fitted by the relation deduced from a 
curling reversal mode in an infinite cylinder |5(| , (3^] : 

h sw (9) = - a(a + 1) (34) 
v /a 2 + (2a+l)cos 2 (0) 

The single adjustable parameter a = —k (Ro/r) 2 is defined by the geometrical parameter k Q], by the exchange 



length R = 20 nm , pi 



and by the radius of the wire r. The experimental points H sw (9) are fitted in Fig. ||. 




FIG. 2. circle: measured position of the switching field Hsw for different angle of the applied field. Line: One parameter fit 
with the curling formula Eq. (^) 

We obtained a = —0.15 (which corresponds to r of about 60 nm). The relation between the angle of the applied 
field 6 and the angle of the magnetization ip c is: 

tan(9) = ^ltan(ip c ) (35) 



7 



The curve Ah, evaluated from Eq. (|3^ ) by numerical resolution with a polarization in the direction of the wire 
axis e p .v = sin(tpo), is plotted in Fig. |3| and Fig. ^, together with the experimental data. A strong discrepancy from 
the linear curve of Ah(I^) at small current pulses can be observed. Above a critical current corresponding to about 
10 7 A/cm 2 the linear fit gives a parameter c = 190, which is in accordance with the rough evaluation of Eq. Q33) . 
This critical current below which the linear regime failed in Fig. [I| could be interpreted following Ref [Q and Ref || 
as the current needed in order to excite spin waves or other magnetization inhomogeneities (37| . The curve given by 
Eq. (|3^) can then be plotted without adjustable parameter (Fig. The divergence at 90° is due to the numerical 
resolution of Eq. (|3^) (numerator and denominator tend to zero at 9 = 90° angle). 




1(A) 



0.0002 0.0004 0.0006 0.0005 0.001 0.0012 



FIG. 3. Parameter Ah — AH m ax(Ie) / H a as a function of the pulsed current amplitude. jj,oH a = 300mT is the anisotropy 
field. The linear fit (continuous line) gives c w 190 (see Eq. (^2|)). The dashed line is the maximum magnetic field induced by 
the pulsed current 




VI. CONCLUSION 



A systematic thermokinetic description of a metallic ferromagnetic layer open to electronic spin polarized reservoirs 
has been performed. At constant temperature, assuming the two current approximation and neglecting direct action 
of the magnetic field on charge carriers, five coupled transport equations account for the complexity of the system. 
The approximation of the explicit uncoupling of the transport processes leads to the known results about GMR and 
Landau-Lifshitz-Gilbert equations for magnetization dynamics. Within this approximation and on the basis of the 
conservation equation of the magnetic moment, the description of both polarized current and magnetization dynamics 
leads to a generalized Landau-Lifshitz-Gilbert equation. The application of this model to experimental data about 
current-induced magnetization reversal is performed. The existence of a critical current indicates that the kinetics 
of magnetization inhomogeneity also plays an important role. However, the comparison with experimental data 
shows that the derived thermokinetic generalized Landau-Lifshitz-Gilbert equation provides a description of the basic 
mechanism responsible for the effect of polarized-current-induced-magnetization reversal. 
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VIII. APPENDIX 



The appendix is structured in three parts. In the first part (A), the equation of for the difference of chemical 
potentials A = A[i is derived from equation ( |l7j ) in the case of steady states. In part (B) the equation is applied 
to the simplified case of GMR or spin accumulation, where the polarization axis is assumed constant through the 
interface. The GMR of the interface is deduced. In part (C), the equation is applied in the framework of the 
experimental study of polarized-current-induced-magnetization-reversal, where an abrupt charge of the polarization 
axis occurs at the interface. 

(A) Assuming that the kinetic coefficients which coupled the dynamics of the magnetization and the electric currents 
vanish, we obtained the following set of kinetic equations (|lj]) : 



Jn 




Jp 






e 



1 -0 

-(3 1 




-^L c 

(T C 



' dz 



A 



(36) 



In the stationary state -^Jn(z) = 0, and assuming that (3, ctq and 7 are approximately independent of z the 
diffusion equation of the chemical affinity is deduced 



d 2 Ho 







d 2 A 



d(jA) \ 
dz J 



(37) 



Inserting in (|3q) yields 



dJr, 
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d 2 A 

dz 2 



0(7 A) \ 
dz J 
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and by integration, 



Jn 



(TO 
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(P 2 
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dA 

dz 
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(39) 



where we assumed that J p (oo) = . 

On the other hand, from the conservation equations we have 



d 



;(Nl 



N ) = I 



k — l — >/c jk- 

V 



■ k+1 



2<i> fc 



2 7 (fc-l,fc)J p 



(40) 



At the continuum limit, we obtain the following relation 



drip 
~dt~ 



2L CC A - 2 7 J p 



(41) 



where n p is the density of spin polarized conduction electrons. 
Equation @ rewrites 



dJp 
dz 



2L CC A - 2 7 J p 



drip 
~dt 



(42) 



where -S^ is constant for steady states . Furthermore, inside the ferromagnet and far away from the interface, J p 
is constant, whence 



dn p 
~dt 



9r°° A 







(43) 
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where we assumed for simplicity that U^ c = 0. Together with (|42|), ( p9[ ) and (B8|) the differential equation for A(z) 
is obtained: 

where the spin diffusion length l s f is given by 



. Mi-/? 2 ) 

= \ 2eL M ' (45) 



the "Domain Wall" diffusion length J^vc is given by 



low - J^P, (46) 



and the parameter k is given by 



2B 2 

k s ijrws (47> 

(B) APPLICATION TO OMR. 

Assuming that the polarization axis is the same for all sub-layers, we have low — and the last term in ( |i^ ) 
vanishes. 

The chemical affinity obeys the diffusion equation: 

d 2 A 1 

The chemical affinity and the total chemical potential are then 

A(z) = ae 1 ^ + fce^f (49) 



/Lto(z) = d + cz+ /3A(z) (50) 

where a, b, c and d are constants. The electric field E(z) is defined by —eE(z) = SM. so that c = ~eE{oo) — 
—eJpf/aQ. Under the condition of continuity of the currents of the two spin channels at the interface (no surface 
scattering) : J±(0~) = J±(0 + ) we have a = ^f^r\ Jn- The spin polarized current on the left side of a single 
interface (b=0) is deduced: 



Jp{*) = ^(%-P^) =JN0[e-*-l) (51) 



The electric field fgf = -E(z) is: 



• \ dz dz 

duo — 



and the supplementary potential due to the spin-polarized current is 

eJw\ 13 

dz = 2- ..„, 

cr / 1 - /? 2 cro 



AV = I + °° f - ^ = 2-^-, (53) 



from which the GMR resistance ( j22| ) is deduced. 
(C) APPLICATION TO SPIN TRANSFER 

In the case of an interface composed by an incident current of conductivity asymmetry (3, polarized in the direction 
e p , entering in a ferromagnetic layer polarized in the direction u with conductivity asymmetry —f3, the Equation (|49|) 
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and (^0() still hold in the left and right hand sides of the interface. However, the change of the polarization axis at 
the interface leads to modify the continuity equation of the current of the two spin channels J±(0~) = J±(0 + ) =p jJ p - 
The integration constant now reads 

<■ - eW/3(1 7 } Jn (54) 



a (l-/? 2 ) 

and the expression of the polarized current is 



J p (z) = J N (3{1 - 7) I e W - l ) (55) 



a ' e-mail : j ean-eric .wegrowe® ip e . dp . epfl . ch 
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